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ABSTRACT 

Whenever  active  or  passive  devices  are  irrcluded:  irr  air 
electronic  circuit  that  is  to  be  analyzed:  by  a;  ciimputer , 
appropriate  models  for  these  devices  imxs?t  bee  djeewedJocpesdc. . 
The  majority  of  the  models  that  Enave  beerr.  adapted:  fxxrr  airar- 
lyzing  a  circuit  on  a  computer  have  negrLecrted:  the:  edzferr:  of: 
ambient  temperature  on  the  element.   Temperature  var±atierc 
has  been  shown  to  have  hazardous  effects  an  electroaiic: 
components.   This  paper  deals  writh.  the  medif^catrLxoi.  o£"  semee 
existing  models  to  include  temperature  effects^  ani  s&xxwss  the: 
results  of  such  modifications.   It  is  done  far  six.  basic: 
electronic  components  (i.e.,  resistor,  capacitor,  inductor, 
diode,  bipolar  transistor  and  field-effect4  transistor) . 
Suggestions  for  extension  of  this  work  are  also  included. 
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I.       EffEHODUCTTON 

A.  GENERAL 

It    is    necessary   far   the  electrarrix:  cbessiigjr.  eaigdneerr  to: 
predict   and   identify   tlie  respxnr^e  a£  el  eclr-unix:  sysitemss  to: 
temperature   effects   caused  by  the:  dezviTrer  iteeeLr:  orr  thef 
environment.      Therefore,,    in.  the  s?tndy   off  sejnixxmrductorr 
electronics    the   desirability   of  computer:  sdmuLatlim- becomes 
readily   apparent.      The   equations^  desscribirrg:  semiconductor: 
behavior   in   all   states   are   complex   ami  nojilineax:; :  hence^ . 
solution   by  manual   means    is  diiFffcuXt:  even  w hear  many 
simplifying   assumptions    are   made.. 

B.  BACKGROUND 

1 .   Circuit  Analysis  by  Computer 

When  a  computer  is  used  to  analyze  an  electronic 
circuit,  a  sufficient  amount  of  information  concerning  the 
circuit  parameters  must  be  supplied  and  device  terminal, 
relationships  in  the  network  equations  describing  the 
circuit  must  be  determined.   This  is  not  a  difficult  task 
for  passive  components  such  as  resistors  and  capacitors. 
Active  devices  may  be  described  to  computer  analysis  pro- 
grams in  terms  of  models  which  mathematically  simulate 
the  behavior  of  the  devices..   It  has  become  increasingly 
important  to  develop  accurate  models-  of  these:  devices-  for" 
computer  programs  used  in  design  and  analysis  of  such  circuits 


Analyzing  the  networks  requires  that  thee  dErnaredc 
equations  must  be  in  some  particular  farm  ami  tfaedirr  paramee-- 
ters  specified  properly  by  the  user.   The?  c±rcu±t  edJemeirts; 
used  in  a  network  could  include  lumped  resistances,,  induc- 
tances, and  capacitances,  current  and  voltage  saouross? 
(dependent  and  independent)  r   and  active  device?  mmfelsf  (i:.e^., . 
diodes  and  transistors) . 

In  the  application  of  the  Computer  Aided  Analysis: 
Programs  the  need  often  arises  for  a  circuit  element,  that: 
is  not  included  in  the  list  of  standard  elements..  Thee 
necessity  of  using  linear  hut  nonstandard,  elements;  (sitchr. 
as  small- signal  transistors  or  ideal  transformers)  arisesc 
in  many  programs.   In  the  case  of  linear  nonstandard  ele- 
ments, general  equivalent  circuits  have  often  been  developed 
by  others;  they  can  be  applied  to  programs. 

2.   Simulation  of  Circuits  with  Active  Elements 

Active  networks  contain  dependent  current  and 
voltage  sources  or  negative  resistors.   When  we  attempt  to 
analyze  an  active  circuit  by  computer  we  encounter  some 
problems.   The  first  one  is  to  find  a  set  of  mathematical 
equations  to  relate  the  current  and  voltage  characreri sties 
of  each  device.   These  equations  can  be  obtained  by  one  of 
two  ways ;  they  can  be  derived  from  mathematical  description 
of  the  physical  phenomena  occurring  during  the  operation 
of  the  device  or  might  be  obtained  by  empirical,  observations 
of  the  behavior  of  actual  device ..   The  degree  of  the 
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prediction,    over  a  wide  ramge:  ccf:  opeaiatrixcrn  aitmrtt  axttrail 
behavior   of   the   device,    determines;  ther  vaaHuee  off  tiresse 
equations.      Secondly   the  parameters;  in.  ttressec  equations;  must 
be   known   or   measurable   or   the   required   inxo  rmari-oir  should 
be   derivable    from   terminal  meas^ireiiiejTts..     Thee  octherr  factor: 
is    the    tolerance   of   th.e   parameters .. 


II.       DEVICE -MODELLING"  METHOTJS- 

A.   INTRODUCTION 

The   objective   of  modelLing:   is   tra  pxaxsssi:  frixnr  tiret 
fundamentals,    that    is,    front  martenrial.  pxapeartdjess  andimor^- 
phology   of   a   distributed  device?  tec  a=  nradiedl  ass  pcrese-nTted; 
below: 

(1)  To   obtain   systematically  a   lumped,    network— like 
model    for    a   distributed  device,,   a)    one  sinxuLcL  writer 
equations   which   are   distributed  airaXugs~  roc  oenree  off 
Kirchhoff's    two   network    laws,,   the  currentt-Tioda  law 
or   the   voltage- loop    law.      Since    the   desired  modeil 
is   network- like ,    it   must   obey  the    two    lajvs    or 
Kirchhoff.      b)    One    should   remove,    by   some   method 

of  approximation,  all  aspects  of  the  mathematical 
description  which  can  be  traced  to  the  distributed 
nature  of  the  device. 

(2)  To  obtain  models  which  describe  the  transient 
performance  of  a  device,  one  might  take  advantage  of, 
if  possible,  knowledge  of  its  steady-state  properties. 
The  reason  for  this  is,  of  course,  that  many  aspects 
of  steady-state  behavior  are  much  easier  to  deter- 
mine than  their  transient  counterparts. 

In  practice,  a  model  is  rareXy  perfect..   Thus;  we 
should  search  for  a  model  whose  elements  provide  the  best 
combination  of  linearity  and  frequency  independence. 
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In   this    section   the   four  foasic  ag^nccacheEss  will  hse 
presented   for   developing   equivalent   circuits  for  esLemerrrs: 
that   are   not    standard   in   computer  programs.      Thees^ee  arrproa:che5: 
can  be   used   individually   or   can   be   combined   to   farm  more 
complicated   device  model. 

B.  PIECEWISE  LINEAR  MODELLING 

In  piecewise  linear  modelling  technique  nonlinear 
elements  can  be  analyzed  by  replacing  all  the  nonlinear 
components  with  equivalent  circuits  r   constructed,  with 
the  linear  elements  recognized  by  the  computer  program.. 
The  nonlinear  characteristics  are  averaged:  over  tire  range: 
of  interest  and  represented  approximately  by  linearized 
characteristics.   For  devices  passing  from  one  operating 
region  to  another  (such  as  a  transistor  moving  from  the 
active  region  to  saturation),  a  different,  linearized 
equivalent  circuit  can  be  proposed  for  each  region. 

C.  MATHEMATICAL  MODELLING 

Those  circuit  elements  that  are  not  standard  (defined 
as  a  circuit  element  that  is  not  included  in  the  list  of 
built-in  models  of  any  computer  library)  can  be  modelled 
using  only  the  mathematical  relationship  which  describes 
their  terminal  characteristics.  The  procedure  to  follow 
in  developing  equivalent  circuits  using  mathematical  model- 
ling technique  is  quite  straight  forwards  In  many  cases-,, 
the  modelling  problem  can  be  reduced  to  the  synthesis  of 
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a   simple    RLC  network   sometimes   rexpxiirxirg:  three  usaer-  off 
dependent    current   sources.      In  otlieir  cases?,,  3L  Hiitttfl3se 
inventiveness    is    useful. 

The   equations    describing   the    components   must  bo  Linear 
since   only    linear   elements    are   used  in.  three  cfev.eianmen t:  of: 
the   model.      The   availability  of  R.,.   L„   ami:  C  edJeanerrlrsf  in: 
the   construction   of  the   equivalent  circuit:  ppeaandittss  tiree 
modelling    of   linear   differential   equations  with    consitant: 
coefficients . 

D.      PHYSICAL   ANALOGY:      PIECEWTSE   LINEAR  MODELLING 

The   modelling   of   a  nonlinear  element:  with:  air.  algobxaic: 
voltage-current    relationship    (sacii  as  nonlinear  resistance) 
was    discussed   in   Section   B.       In   this   section    the  techniques 
for  modelling    a   nonlinear   reactive    element   with    a  non- 
algebraic   voltage-current    relationship    is    discussed    (such 
as    inductance    or    capacitance).      For   better   understanding 

* 

one  may  refer  to  examples  given  in  Ref.  17.   Since  non- 
linear inductance  and  capacitance  are  algebraic  functions 
of  magnetic  flux  and  charge,  respectively,  which  are  not 
recognized  by  the  program,  a  reasonable  solution  to  this 
modelling  problem  is  to  find  a  physical  analogy  which  will 
allow  flux  or  charge  to  be  treated  in  those  units  recog- 
nized by  the  program. 

The  physical  analogy  acts  as  a  black  box  which  converts 
the  physical  relationship  into  terms  which  are;  accopirable 
to  the  program  as  shown  in  following  algorithm.   Therefore, 
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elements    such    as   nonlinear   imdcicttarLces  and: 


;s  cair. 


PHYSICAL 
COMPONENT 

—9 

PHYSICAL 

ANALOGY 

-9 

COMPONENT  EQUIVALENT 
CIRCUIT  USING 
V-I  RELATIONSHIP 

— 9 

COMPUTED 

• 

PROGRAM 

Figure    2-1.      Physical   Analogy  Concept:  Usae 
in   Physical    Component  ModelLing:.. 

be  modelled   if   an   analogy  between   the  cdiairact^r±S5t±c5i:  to:  ber 
modelled   and   a  voltage-current   relationship   acceptabLe    to 
program   can   be    established. 

E.   PHYSICAL  ANALOGY:   MATHEMAT I CAL   MODELLING: 

Almost  all  of  the  Computer  Aided  Circuit  Analysis: 
Programs  perform  calculations  in  terms  of  voltages  and 
currents  only.   The  parameters,  such  as  charge  and  magnetic 
flux,  have  not  been  included  in  the  programs  and  cannot  be 
used  directly.   If  such  conditions  exist,  it  is  necessary 
to  develop  methods  for  expressing  the  new  parameters  (charge, 
magnetic  flux,  etc.)  as  functions  of  voltage  and  current  or 
to  redefine  some  new  elements  in  addition  to  the  elements 
included  in  the  standard  list  that  will  perform  the  required 
analysis  in  terms  of  voltage  and  current. 

The  solution  of  mechanical  and  thermal  problems  by 
electrical  analogy  techniques  is  an  example  of  this  kind. 
Establishing  the  analogy  is  the  first  step  in  the  deriva- 
tion of  an  equivalent  circuit.   The  second  step  is."  Linking: 
the  derived  analog  network  to  the  remainder  of"  the  model' 
which  includes  the  element  terminals.   The  element  terminal 
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characteristic  must   always   be  defined;  in  tterans-  ccfr  viadHta^ec 
and   current.      The    charge-control    transistor-  nratteL  iss  a= 
good   example   of   physical    analogy   techniques   in.  nrodeilXijTg; 
the   other   charge   equations   which   govern   the   transistor- 
response.      See    Ref.    17,    Appendix  A. 2. 
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III.   TEMPERATURE  EFFECTS 

A.   TEMPERATURE  EFFECTS  ON  RESISTOR 

The  resistivity  p  of  a  conductor  is  the  ratio  of  the 
electric  intensity  E  to  the  current  per  unit  cross-sectional 


area. 


P  -  E/(I/A) 


3.1 


The  resistivity  of  all  conducting  material  is  affected  by 
its  temperature.   A  plot  of  resistivity  versus  temperature, 
for  a  metallic  conductor,  is  given  in  Fig.  3.1.   The  curve 
may  be  satisfactorily  represented  by  an  equation  of  the  form 


2      3 
P  =  P0  +  g-,1  +  g9T   +  g^T 


I1  '  62*   '  63 


3.2 


where  pQ  is  the  resistivity  at  0°  C,  g, ,  g?,  ...,  are 
constants  depend  on  the  characteristic  of  a  particular  material 


Figure  3-1.   Resistivity  versus  Temperature 


15 


For   temperatures  which   are  merit  trcrcr  grreartt  tthee  ttBamss  iir. 
T     and  higher  powers   may  be  neglected  aird:  carr.  tree  wrditterr.  ass 

P   =    P0    +   &iT    •  25,-3c 

The    temperature   coefficient   definition  of  resrs^trLviity    r_s; 

a   =   g2/P0    =    (p-pQ)/(P0T)  3E..4- 

Since   the    resistance   of   a   given   conductor/  is  p rupm:tijocrraiL 
to    its    resistivity    the    following   equation   can.  he  wriitten 

R   =    R0(l    +   aT)  33..5.: 

where  R„  is  the  resistance  at  0a  C  and  R  ils^  the  resistance. 
at  T°  C. 

The  resistivities  of  nonmetals  decrease  with  increasing 
temperature,  and  their  a's  are  negative  Electrolytes  also 
have  negative  a's. 

B.   TEMPERATURE  EFFECTS  ON  INDUCTOR 

Temperature  variation  is  usually  accompanied  by  a  change 
in  inductance,  resistance  and  self-capacitance  of  a  coil. 
The  last  two  types  of  change  are  not  important  in  compari- 
son with  the  change  in  inductance.   The  formula  for*  single- 
layer  solenoid  given  by  Sturley  [Ref.  6]  is 

1  =  (r2N2)/(9r  +  1(H)  3.6 

where  I  =  length  of  winding 
r  =  radius  of  winding 
N  =  total  turns  in  the  coil- 
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The   change   of  the   induxrtamre?  fair  ar.  iircneaasaee  iir,  temperature: 
of  T°    C   can  be  shown  to  free   [Ref..   66]" 

AX   =  bL(AT)  3:7" 

where  AL=    change    in   inductance 

b   =   coefficient  af  Lirneasir  expmrsdxnT.  off  thee  conductor: 
AT  =   change  in  teanp:en?atnre£.. 
Equation    3.7    is   valid:  when,  the;  oooefff  i  ci  errtsf  off  rani  ail  and; 
axial    expansion,   b      and  br,.  are:  approximately   equals  . 

It    is   possible    to  make   the  inductance   L  independent; of 
temperature  by   satisfying:  Eq..   35..8i,.  the:  ratio;  of;  axial!  to; 
radial   expansion. 

hff/b      =   Z  +■    (9r/l(U)  3^8 

By  choosing  the  coil-former  and  the  conductor  material 

4 

properly,  the  temperature  coefficient  for  inductors  can  be 
minimized. 

C.   TEMPERATURE  EFFECTS  ON  CAPACITOR 

Temperature  variation  has  two  effects  on  capacitors, 
namely  to  change  the  dielectric  permittivity  and  the  dispo- 
sition of  the  plates.   Some  typical  capacitance- temperature 
coefficients  will  be  given  in  Chapter  IV.   Permittivity  of 
the  insulating  material  is  considerably  changed  with  tempera- 
ture.  Capacitances  for  parallel  plate  and  cylindrical  types 
are  given  with  the  following  equations, 

parallel  plate::  C  =  (A..e;)/ci  319. 


IT 


cylindrical:   C  =  CZtit.  JL.E))/Tccg(:av'Ir)         SIDC 

where  area  of  plate  (A),  and  Length,  of:  canac±tDcr  (&;),,,  are-: 
temperature  dependent  also,  but  this  dependence  can  be 
negligible  when  compared  with,  the  dependence  off  p;ermi3rtivity , 
In  both  cases  capacitances  are  "the  function,  crff  pxa^miittivity 
e.   In  addition  to  temperature  deprendencer  afr  c,.  it  is? 
dependent  of  frequency,  but  tie  latter  efffeeczt  is?  a?,  ddiffeorent: 
story. 

Improvement  of  the  manufacturing;  techniques-  minimizies- 
the  temperature  coefficient  for  capacitors.. 

D.  TEMPERATURE  EFFECTS  ON  DIODE 
1.   Introduction 

There  are  two  kinds  of  charge-carryings  particles"; 
one  is  negative  (the  free  electron)  of  mobility  u  ,  and 
the  other  is  positive  (the  hole)  of  mobility  \i    .   These 
particles  move  in  opposite  directions  in  an  electric  field 

E,  but  since  they  are  of  opposite  sign,,  the  currents  are  in 
the  same  direction.   Hence  the  current  density  J  is  given  by 

J  =  ^nun  +  pup)  G  E 
where  n  =  magnitude  of  free-electron  concentration 

p  =  magnitude  of  hole  concentration 

e  =  electron  charge  in  coulombs. 
For  an  intrinsic  semiconductor  n  =  p  =  n._,  where  n.  is  the 
intrinsic  carrier  concentration-   The  intrinsic:  carrier: 
concentration  n-  is  temperature  dependent,  such  as 
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n?  =  AT3exp(-EG0/kT)  . 
It  has  been  shown  experimentally  that  E~  depends  upon 
temperature  also  as  pointed  out  below 


EG  =  EG0  "  qT 


where  Ern  is  the  energy  gap  of  the  semiconductor  at  0°  K 


'GO 


-4 


and  q  is  a  constant  (q  =  2.4  x  10    eV/°K  for  silicon  and 
q  =  3.9  x  10    eV/°K  for  germanium).*   Such  a  large  change 
in  conductivity  with  temperature  places  a  limitation  upon 
the  use  of  semiconductor  devices  in  some  circuits. 
2.   Temperature  Dependence  of  Parameters 

For  a  p-n  junction  the  current  I  is  related  to  the 
voltage  V  by  the  equation 


I    =    Io[exp(V/mVT)-l] 


3.11 


Temperature  dependent  parameters  of  the  Eq.  3.11  are  I 


and  V~T. 


I  =  -I 


n 


V 
(a) 


Voltage 


2  5°C 


(b) 


1  Micro  amp. 

- @  25°C 
100  micro  am1 
@  150°C 


Figure  3-2.   Typical  Reverse  Biased  Diode  (a), 
and  its  Characteristics  (b) 


:Strictly  speaking,  q  is  not  a  constant.   For  example,  fo: 
silicon  q  =  2.4  x  10    eV/°K  only  when  T  >  200°  K.   For 
germanium,  q  =  3.9  x  10"4  eV/°K  only  when  T  >  150°  K. 
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I      is    called   the   reverse   satunravfcicm!  (ainnraant!  andi  gdiveir  by 
the   equation 

I      =    K  Tn   exp(-VGQ/mVT)  37.  IX.  a 

for  Germanium:   m  =  1    n  =  2    Y^  =  Q..78l5r  vcrilts- 

for  Silicon:     m  =  2    m  —  1.5  v\~-  =  L..2IL  volts 

From  experimental  data  we  find  that  th-e?  reverse  ssasturatlonr 

current  increases  approximately  seven  per  cent:  pear  °:  C 

Since  (1.07)    =  2.0  reverse  saturation  current:  approximate:-- 

ly  doubles  for  every  10°  C  rise  in.  temperature..  V.-^  iss  thee 

"volt  equivalent  of  temperature"  defined:  by 

VT  =  T/II600  . 

Another  important  parameter  for  diode  ts   diffusion  or:  storage 
capacitance  and  given  by  * 

CD  =  t  I/m  VT  1.12 

where  t  is  the  diode  time  constant  equals  the  mean  lifetime 
of  minority  carriers. 

Reverse  leakage  resistance  also  has  some  temperature  depen- 
dence; since  it  can  change  depending  on  the  manufacturing  pro- 
cess, no  generalized  temperature  coefficient  can  be  stated. 
It  should  be  determined  for  each  family  of  the  diode. 

The  temperature  coefficient  of  the  diode  parameters  can 
be  derived  in  differential  form  from  Eqs .  I. 11  and  3. II. a 
by  taking  the  derivatives  and  making  some  approximations 
to  yield  the  following  equations: 
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(l/I0)(dI0/dT)    =    foflty  *   0%,/Wroy  3313: 

dV/dT  =    [V-CVG0   *■  rairVr)l/T  33114- 

at   fixed   current. 

E.      TEMPERATURE  EFFECTS  <M  TRANSTSTQR 
1.      Introduction 

A  junction   transistor  consists-   off  a,  sdiDxixnr.  orr 
germanium   crystal,    in  which  a  Tayerr  of  n.-typa  orr  p7-type: 
material    sandwiched  between  two   iayerrs  off  p:— type:  orrnrtype: 
material,    is    called   a  p-n-p   cnr  n-p/-n  trnarrsistDrr. .     Fbrr  ai 
qualitative   understanding   of  the   faniL  af"  the   input  and:  out"-- 
put    characteristics,    one    can   consider   that   the    transistor 
consists    of    two   diodes    placed   in   series    "b<ack    to   back."      It 
has    four   operating   regions    defined   as    follows: 

Region  (a) :   Cut-off  region  (emitter  reverse  biased, 

collector  reverse  biased) , 

Region  (b) :   Normal  transistor  operation  (emitter 

forward  biased,  collector  reverse  biased) , 

Region  (c) :   Saturation  region  (emitter  forward  biased, 

collector  forward  biased) , 

Region  (d) :   Inverse  transistor  operation  (emitter 

reverse  biased,  collector  forward  biased).   The 

input-output    characteristics   of  the    transistor  in 

different   modes    are   given   in  many   texts. 
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2.      Temperature   Dependence  of  Faxametexs: 

All   of   the  parameters  of  ijr.tjerresit  in:  c^mp-uterr 
application   concerning    transistors  change:  withr.  temperature * 
Some   of   them   are   very   sensitive    to   temperature   changes-, 
while   others   may  be   assumed  essentially-  cnnsftairtt  inn  the: 
analysis.      The   most   tentperatuxe  SECTraErLtrfiv.ee  paxameteTSs  aree 
saturation   currents,    Ip^   and   Ip^,.   cumrent;  ga±rr  factors 5, . 
aN   and   aT,    and    leakage   resistances,.  Kp  and-  RT^„.     Parameters: 
like    time    constants,    t^,    z-rr   cap aciitancjes;  Cp-,.  Cp-j.  pxopor^-- 
tionalities,    mp,   mr,    change   only   slightly   overr  the: 
temperature   region   of  interest  as  high,  ass  L23;°:  d.. 

(a)    Reverse   saturation  currents:      The  emitter: 
junction   diode   saturation   current    Ip~   and   collectorr  junction 
diode   saturation   current    1^-    are   the   most  .temperature    sensi- 
tive  parameters.      Suppose    that    the    reverse   saturation 
current    at    room   temperature    (300°    K)    is    given   by   Eq.    3. 11. a; 
it   can   be   written   at    temperature    T   +    AT,    as    follows: 

IC0   =    Kl   Tn   exP(-EG0/kT) 
IC0T   =    KXCT   +    AT)n   exp[-EGQ/k(T   +    AT)]  35.15 

equivalently 

I-E0T  =  K2(T  +  AT)n  exp[-EG0/k(T  +  AT)]        3.16 

Assuming  AT  <<  T  and  expanding  the  Eq..  5.  LS:  in  Taylor  series, 
the  following  can  be  obtained: 
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XCOT   =    !C0{1    +KEGO/kTJ    *  mI    W^B    "  %dED1  -k  WJfflM]     3al7' 
After   numerical    substitution   it   tuiirns  out:  tar  lie; 

ICOT   S    IC0(^1   +   °-157   AT) 

where  X  =  0.157/°C  and  is  defined  as  the  treairrreErra^turee 
coefficient  for  silicon  transistor.   Eonr  grrealierr  accuracy , . 
over  larger  temperature  ranges r    the  eaTia..t±orr.  rjiczLuding: 
the  exponential  factor  directly  should  be  usect. 

(b)  Current  gain  factors:   The  current  ga±rr.  factors; 
in  normal  and  inverse  modes  for  common:  base  cxm£rguratxoar 
are  given  below: 

aN  =  hFE/(1  +  hFE}   2nd   aI  =  hFEr/(L  *  hW 
where  hpF  is  normal  mode,  common  emitter  curreirtL  gain,,  and 
hpFT  is  inverse  mode,  common  emitter  current  gain.   otN  and 
aT  do  not  change  appreciably  with  temperature  but  (1-oc^)  or 
(1-a-r)  plays  an  important  role  in  the  following  formulas: 

hFE  =  aN/^1'CL^      and  hFEI  =  ai/(l-ai^     • 
For  example,  the  variation  of  common-emitter  current  gain 
hpp,  with  temperature  normalized  to  unity  at  T  -  25°  C 
for  the  2N524  transistor  is  shown  in  Fig.  3.3. 

(c)  Reverse  leakage  resistance:   Reverse  leakage 
resistance  is  temperature  dependent  to  some  extent.   No 
generalized  temperature  coefficient  can  be  stated  since  it 
varies,  depending  upon  the  process  of  manufacturing..  There -- 
fore,  it  should  be  determined  for  each  family  of  transistor. 
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F.   TEMPERATURE  EFFECTS  ON  FIELD-EFFECT  TRANSISTOR 

1 .   Introduction 

Field-effect  transistor  (FET)  has  two  primary 

operating  regions,  namely  the  triode  region  in  which  the 

device  behaves  as  a  voltage-variable  resistor,  and  the 

pinch-off  region  where  the  current  is  determined  mainly  by 

the  gate  voltage.   The  two  regions  are  separated  by  the 

locus    of   pinch-off   points,    i.e.,    the   points   where   V~DC.    = 

V  -  V~c.   V^c ,  Vn  and  V~c  are  defined  drain  to  source 
p     Go     Do 7       P       bb 

voltage,  pinch-off  voltage  and  gate  to  source  voltage, 
respectively. 

The  operation  of  FET  depends  upon  the  flow  of  majority 
carriers  only.   Therefore,  it  is  a  unipolar  device.   Since 
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the   majority- carrier  current  decreases  wiithi  t&anpBarsiiirEe, . 
the    troublesome   phenomenon    [Rei:..   #1   is  iratr  enranirrteread: 
with    field-effect   transistors.      This  property   procv.ijdesf 
thermal    stability   for   FET. 

2.   Temperature  Dependence  of"  Parameters 

Drain  saturation  current,.  F^crcr  which",  dsperrds?  upon 
channel  geometry,  impurity  density-  and  charge:  carrrerr 
mobility,  is  defined  to  be  the  channel  current:,.  Inc^ 
at  Vpc  =  -Vp.   Although  I^cc  is  rart  greatly  dependent  on 
V^- ,  it  depends  an  temperature  heavily..  TTeraperatiire:  can-. 
influence  the  drain  current  in  two  ways?..   Ffiinssfcly ,.  aa  changee 
in  mobility  causes  a  change  in  resistance  off  channels. 
Secondly,  a  change  in  temperature  causes  a.  change  in  surface 
potential.   The  change  in  surface  potential  works  in-  opposi- 
tion to  the  change  in  mobility.   Thus,  change  in  surface 
potential  constitutes  a  built-in  gate-voltage  compensation 
of  the  mobility  change  in  the  semiconductor.   The  effect  is 
dependent  on  the  D.C.  gate  bias  so  that  at  a  particular  bias 
the  compensation  may  be  optimized. 

Another  quantity  of  interest  is  the  transconductance. 
The  transconductance,  g  ,  has  the  same  temperature  variation  as 
I   does.  The  principal  reason  for  the  negative  temperature 
coefficient  of  Inc,  is  that  the  mobility  decreases  with 
increasing  temperature. 

As  pointed  out  before,,  the  major  effects  of-.'  tempera-- 
ture  on  a  FET  or  alteration  of  the  average  channel  conductivity 
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Normalized  I^oc:  versus  Tempera- 
ture, at  25°  C 

and  the  gate-channel  junction  barrier  potential.   Since 
conductivity  decreases  with  increasing  temperature,  the 

4 

drain  current  decreases  as  the  temperature  rises.   On  the 
other  hand,  barrier  potential  decreases  with  increasing 
temperature,  causing  the  magnitude  of  the  total  gate-source 
potential  to  be  reduced,  hence  increasing  the  drain  current 
By  suitably  choosing  the  operating  point,  it  is  possible  to 
make  these  effects  canceled,  producing  what  is  essentially 
a  zero  temperature  coefficient  operating  point. 
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IV.   COMPUTER  MODELS 

A.   INTRODUCTION 

Almost  all  the  computer  analysis  programs  are:  cfessigneji 
to  solve  the  linear  and  nonlinear  simultaneous  e^pxa^txansE 
which  characterize  the  mathematical  models  used!  toe  prr^aixct: 
circuit  response  for  an  electrical/temperature  eirviiranrnent:.. 
The  user  codes  the  problem  by  entering  the  circuit  tapalogy 
and  parameters.   Programs  contain  a  standard  (buiXt-in.) 
component  model  library.   Non-standard  models  may-  he  enrr- 
structed  by  using  standard  models.   FORTRAN  algrebraxc 
equations  and  control  statements  may  be  coded  to  modify 
standard  model  parameters,  augment  the  standard  models, 
create  nonstandard  models,  create  additional  time  functions, 
temperature  functions,  etc.   The  circuit  is  entered  into 
the  computer  program  by  using  the  program  standard  part 
models  and  specifying  the  circuit  topology. 

In  this  paper  by  "conventional  (or  non-standard)  model," 
it  means  the  model  of  a  device  which  is  isolated  from  all 
environmental  effects,  and  by  "modified  model,"  it  means  a 
model  with  temperature  effects  taken  into  account.   In 
addition  to  this  task,  some  additions  or  some  simplifications 
will  be  done  when  necessary,  depending  upon  device  characteris 
tics . 

In  order  that  electronic  components  may  be  handled 
quantitatively,  a  set  of  mathematical  relationships  must.be 
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derived  relating  their  terminal  characteristics.   Deriva- 
tions will  not  be  attempted  in  this  paper  since  they  are 
documented  elsewhere  [Refs.  2,  3,  4,  7,  18,  22];  only  the 
results  will  be  presented. 

B.   RESISTOR  MODEL 

1.   Conventional  Model 

The  resistor  model  is  defined  by  the  following 
current  equation 

iR  -  VR/R 
where  VR  is  the  instantaneous  voltage  across  the  resistor  R, 


->  iV 


T,  • AV o   T9 

R  z 

4 

Figure  4-1.   Resistor  Element 


Since  the  element  is  bilateral,  the  choice  of  which  is  con- 
sidered to  be  terminal  1  and  which  is  terminal  2  depends  on 
the  direction  of  current  flow  (from  T,  to  T~) .   This  element 
has  just  one  parameter:   the  resistance  value  R. 
2.   Modified  Model 

Temperature  dependence  of  resistor  was  discussed  in 
Chapter  III,  but  the  Eq.  3.5  can  be  rewritten  in  slightly 
different  form  such  as 

RT  =  R  ±  AR  4.1 

where  AR  is  the  amount  of  resistance  change  due  to  tempera- 
ture change  and  equivalent  to 
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AR  =  ±aR(AT)  4.2 

The  sign  depends  on  material  properties.   Assuming  resistor 
is  linear,  the  following  modified  resistor  model  for 
computer  can  be  used.   The  determination  of  AR  requires  the 
parameter  a,  which  can  be  found  in  many  manuals  or  reference 
books . 


■>1R 


R  AR 


i  • ;  •  '/ '     '  T2 


Figure  4-2.   Modified  Resistor  Model 


C.   CAPACITOR  MODEL 

1 .   Conventional  Model 

A  capacitor,  shown  in  Fig.  4.3,  is  a  "bilateral 
element  with  terminals  T,  and  T~,  having  a  capacitance  C 
The  current  flowing  through  a  capacitor  is  given  by 

ic  =  C[d(V1  -  v2)/dt] 


>xc 


Figure  4-3.   Capacitor  Element 


V IS ^2 
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2.   Modified  Model 

The  more  realistic  capacitor  has  a  shunt  leakage 
resistance  and  a  series  resistance  as  shown  below  that  are 
identified  with  R,  and  R~  ,  respectively. 


o  T. 


Figure  4-4.   Capacitor  Model  with  Resistive  Components 


Temperature  dependence  of  capacitor  was  discussed  in 
Chapter  III.   Very  small  series  resistance  /in  the  order  of 

0.1  ohms)  and  very  large  shunt  resistance  (in  the  order  of 

9 

10   ohms)  can  be  assumed  constant  with  temperature  change. 

Some  typical  capacitance- temperature  coefficients  are  given 
in  the  following  table  and  many  others  can  be  found  in  the 
references  given  in  the  footnote.* 

Table  4-1.   Capacitance-Temperature  Coefficients 
Insulating  Material     Capacitance -Temperature  Coefficients 
Varnished  Cambric     +2100  parts  in  106  per  degree  Centigrade 
Synthetic  Resin       +1600  parts  in  106  per  degree  Centigrade 
Enamel  +470  parts  in  106  per  degree  Centigrade 

Ceramic  +100  parts  in  106  per  degree  Centigrade 


*(1)  Handbook  of  Chemistry  and  Physics;  (2)  International 
Critical  Tables  Vol.  2,  McGraw-Hill;  (3)  Tables  of  Dielectric 
Materials,  Vols.  I-IV,  MIT,  January,  1953. 
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Defining  a  change  in  capacitance  as  AC  the  modified  capacitor 
model  can  be  obtained  shown  in  Fig.  4-5,  in  which  AC  should 
be  calculated  and  added  to  the  model  with  the  aid  of  tables 
(given  in  the  footnote  on  page  30)  and  the  Eqs .  3.9  and  3.10. 
The  change  in  capacitance  is  very  important  if  the  capacitor 
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Figure  4-5.   Capacitor  Model  with  Temperature  Effect 


is  an  element  of  resonance  circuit  because  a  small  change  in 

4 

capacitance  may  cause  a  large  drift  in  resonance  frequency. 

D.   INDUCTOR  MODEL 

1 .   Conventional  Model 

The  relationship  for  current  is  given  by 


iL  =  (1/L)/  (Vj  -  V2)  dt 


where  L  is  the  value  of  inductance  in  henries  and  the  dif- 
ference, (V,  -  V2) ,  is  the  voltage  across  the  element.   The 
inductor  is  diagrammed  in  Fig.  4-6.   Since  the  inductor  is 
bilateral,  the  choice  of  which  to  be  considered  T,  is  com- 
pletely arbitrary.   This  element  has  two  parameters:   induc- 
tance L,  and  Q  (of  the  coil).   Many  applications  of  inductors 
have  practical  equivalent  circuit  that  also  includes  the 
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Figure  4-6.   Inductor  Element 

quantity  Q--not  just  pure  inductance.   Frequency  dependent 
r  is  given  by 


r  =  Q.27rfL 


4.3 


where  f  is  the  frequency  at  which  the  analysis  is  being 
performed.   The  significance  of  Eq.  4.3  is  best  seen  in  a 
sketch  of  resistance  versus  frequency.   The  range  for  the 
Q  model  is  that  r  becomes  unrealistically  small  as  frequency 


r  '< 


is  typically 
pecified  at 
ome  test 
requency . 


Figure  4-7.   Resistance  r  versus  Frequency  f 

becomes  very  low,  and  at  high  frequency  the  model  would  show 
r  approaching  infinity  which  is  also  unrealistic. 
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2.   Modified  Model 

Temperature  dependence  of  inductance  was  given  in 
Chapter  III  with  Eq.  3.7.   Additionally  D.C.  series  resis- 
tance and  the  shunt  leakage  resistance  should  be  involved 
into  the  model,  which  are  assumed  to  be  independent  of 
temperature.   The  schematic  which  takes  all  the  effects 
explained  above  into  consideration  is  shown  in  Fig.  4.8. 
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Figure  4-8.   Inductor  Model  with  Temperature  Effect 

To  determine  AL  one  should  first  determine  the  parameter 
b  in  Eq .  3.7.   The  value  of  b  can  be  found  in  the  handbooks 
specified  in  the  footnote  in  previous  section.   For  better 
understanding  of  importance  AL,  look  at  the  example  given 
below. 

Example:   For  a  tank  circuit,  f   =  1/(2tt/L~C).   If  Af  is 
the  change  in  oscillation  frequency  due  to  AL  change  in  L, 
the  following  relationship  can  be  obtained  after  proper  sub- 
stitution and  expansion  by  the  Binomial  Theorem  where  only 
the  linear  term  is  retained: 
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The   value   of  b    for  capper  is  E-G6  x:  HOC     '  pcaartss  pear  °=  CC .     TEusf 
the    ratio   change   in   frequency"  due:   to:  thee  uirhxxctkexedc  expansion-. 
of   a   copper   coil    is    8   x  10        praxis  pear  <*  C,.  arr  foxr  ATT  ==  3D°    C 
the   ratio   change   of   frequency    is    Z40    parts    in   10    ..     Frequency 
variations    due   to  coil  expansion,  of  Z44I  and:  214-011  hzz  ait  oscil- 
lation  frequencies   of  I  and  10  Mfcz:,.  r^s^jexxively , .  wouldibec 
expected   under   tliese   conditions.. 

E.       DIODE   MODEL 
1.      General 

Any   diode   can  be  moxfelled:  using:  the;  fodlbwing:  approach . 
The   quantities    of  interest  are:  electxlcall  chaxaxtexi sties: 
which   can  be   observed   at   the   external  terminals-  ( i:.  e:. ., .  a^ 
"black -box"    approach)  .      This   type   of  model. Ting,  requires: 

4 

finding  a  mathematical  expression  to  represent  each  electri- 
cal characteristic.   In  this  paper  this  method  will  be 
emphasized.   Physical  constants  will  be  used  whenever  possible, 
but  if  a  characteristic  cannot  be  related  physically,  an 
empirical  description  will  be  used. 

Based  on  the  operating  range,  diode  models  may  be 
classified  into  two  groups:   "small-signal"  and  "large-signal" 
models.   The  first  type  consists  of  finding  an  operating 
point  and  assuming  that  all  responses  are  linear  about  that 
point.   Therefore,  this  model  is  valid  over  only  a  limited 
operating  range,  namely,  where  the  response  characteristic 
is  linear.   Large-signal  model  is  one:  which  can  be.  used,  over: 
a  wide  range  of  operation.   By  its  nature,  it  is  the  best 
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suited  model  for  use  in  general  osmprutenr  anaXLyssiss..  Tfreenee 
are  three  most  widely  known  models:   (X)  En:erns5-MoXX  [Re£^.  1D1  , 

(2)  Beauf  oy-Sparkes  charge -control  model  [Efefr..  XX]",.  and: 

(3)  Linvill  lumped  model  [Ref.  12],   In  this  paper  Ebers- 
Moll  model  will  be  examined  in  irare  detail  ami  nexEESSsary 
modifications  will  be  given.   The?  Treason:  fmr  th±S5  pr^£er©nr£r 
is  this  model  describes  the  "Id lack -box4''  frenaviarr  with:  redaa-- 
tively  simple  mathematical  equations .. 

2.   Model  without  Temperature  Effects 

The  current  equation  far  diode  was  given,  fredibre:  in 
Eq.  3.11.   Another  diode  current  equation  will  free  oiitnined: 
by  discussing  the  various  factors  in  the  semiconductor:  diodes 
which  have  to  do  with  its  electrical  performance. 

,>         With  the  application  of  an  external  reverse  bias, 
the  p-n  junction  width  will  expand  because  more  donors  and 
acceptors  will  be  uncovered.   This  increase  in  uncovered 
charge  with  applied  voltage  may  be  considered  a  capacitive 
effect,  named  transition-region  capacitance. 

If  bias  in  the  forward  direction,  the  potential 
barrier  at  the  junction  is  lowered  and  holes  from  the  p-side 
enter  into  n-side.   Similarly,  electrons  from  the  n-slde  move 
into  p-side.   A  change  of  voltage  on  the  p-n  junction  will 
alter  the  amount  of  stored  carriers.   This,  then  looks  like 
a  capacitor  to  the  accompanying  circuit,  named  storage  or 
diffusion  capacitance. 

By  taking  the  leakage  resistance  into  account;,  the  follow- 
ing diode  current  equation  can  be  written  as 
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ID  =  I  +  E/RL  +  CD(dE/dt)  +  CT(dE/dt) 


4.4 


These  previously  derived  properties  are  now  all  combined  to 
provide  a  composite  equivalent  circuit,  or  computer  model, 
of  the  p-n  junction  diode,  shown  in  Fig.  4.9. 
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Figure  4-9.   Diode  Computer  Model 


•  Cathode 


The  definition  and  determination  of  the  model  parameters 
will  be  discussed  in  the  following  section. 
3.   Determination  of  Model  Parameters 

Ideal  model  parameters  can  be  read  directly  or 
extracted  from  the  device  specification  sheet.   Unfortunately, 
this  cannot  always  be  done.   Most  data  sheets  do  not  contain 
sufficient  information  to  define  the  entire  model. 
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In  Fig.  4-9,  RR  represents  the  semicomducitxir-  hxLllk;  resiisc— 
tance,  RT  the  junction  leakage  resistance,  CL^  the?  tzran-ssiitniorr 
capacitance,  and  I  is  an  ideal  current  source  which  djescriiiesc 
the  forward  current  as  a  function  of  junction  voltage  E..  Et. 
is  possible  to  write  the  Eq.  3.11  in  slightly  diffTerrerrt:  fornr 
as   given  below: 

1=1    [exp(eE/rakT}-l]  41..5E 

where    e,    k,    and   T   are   electronic    charge  r    Baltzmamr's  constant:, 
and   absolute   junction   temperature    in    aK„    respectively..   L 
and  m   are    constants   which   define   the   shape  of  the  V-E  ciraxac:-- 
teristics.      For   D.C.    operation,    the   network   equatxccrrs;  arre::: 

VD   =    E    +    ID   RB  «L« 

ID   =    I    +    E/RL  4.7 

A  plot  of  the  static  V^  -  Iy.  characteristic  for 

positive  In  is  shown  in  Fig.  4-10.   The  current  is  plotted 

on  a  logarithmic  scale  to  illustrate  that  over  a  wide  range 

in  the  forward  operating  region,  V"n  =  E.   For  determining 

the  unknown  model  parameters  I  ,  m,  and  Rn  one  method  is 

r  O  B 

to  curve- fit  data  points  on  the  forward  diode  characteristic 
to  Eqs .  4.6  and  4.7,  where  1^  is  assumed  to  equal  I.   Another 
approach  to  finding  the  parameters  involves  the  same  assump- 
tion but  not  the  curve  fitting.   One  may  choose  two  points 
on  the  linear  portion  of  the  curve  in  Fig.  4.10  and  label 
them  (I,,  V,)  and  (I-,  V"2).   In  this  range  maybe  expressed 
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with  following 


ID  =  IQ  exp(6VD/m) 


4.8 


where  6  =  e/kT, 


ID  (amps) 


Log 
scale 


10-2H 


10 


10 


10 


X3RB 


ID  =  IQ  ex^  (6VD/m) 


0 


Figure  4-10 


0.5 


1.0 


VD  (volts) 


Diode  Forward  V  -  I  Characteristic  (Typical) 
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Taking  the  logarithm  of  both  sides  and  substituting 
(I,,  V,),  (I2»  V-) ,  then  solving  for  m  gives 


m  =  6(V2  -  V1)/ln(I2/I1) 


4.9 


Substituting  either  data  point  and  the  value  of  m  found  in 
Eq.  4.9  into  Eq.  4.8,  the  following  can  be  obtained: 

I   =  I1/exp(eV1/m)  =  I2/exp(6V2/m)  . 

The  difference  between  the  linear  extrapolation  and  the 
measured  point  is  the  drop  due  to  RR.   Solution  of  Eq.  4.6 


results  in 


RB  "  (V3  "  V/T3 


4.10 


where  E3  =  (m/6)  In  [(Ij  +  I0)/IQ]  £  Cm/8)  In  U3/I0). 


(I4.V4) 


Figure  4-11.   Diode  Reverse 
V  -  I  Characteristic 


From  the  reverse  characteristic  shown  in  Fig.  4-11,  the 
ohmic  leakage  resistance,  RT  can  be  found  to  be 
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RL  =  V^4  +  V  • 


4.11 


There  remain  two  more  parameters  to  be  determined,  namely, 
the  junction  transition  capacitance  CT  and  diffusion  capaci- 
tance C^  in  Fig.  4-9.   A  typical  curve  of  junction  capacitance 
versus  reverse  bias  is  shown  in  Fig.  4-12. 


>  V 


D 


Figure  4-12.   Junction  Transition 
Capacitance  Characteristic 


The  analytical  expression  that  describes  CT  is 


CT  =   K/CVj  -  V1) 


n 


4.12 


where  K,  n,  and  V"T  are  junction  capacitance  proportionality 
constant,  grading  constant  and  junction  contact  potential, 
respectively  (V,  is  approximately  0.9  or  1.0  volts).   For 
VD  =  0.0,  CT  =  C   =  K.   Therefore,  at  any  point  (C,,  V^) 


ci  ■  V(vj  -  V 


Solving  Eq.  4.13  for  n  results  in 

n  =  ln(C0/C1)/ln(VJ 


n 


4.13 


V  • 


4.14 
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Another  way   of   expressing   the   Eq.    5.12  imay  lie  tdte?  fictilll3ro«djffgg: : 

Cn   =    6(1    ♦    I    )/m2Trf  44.JJ55 

u  o 

where  the  parameter  f  is  diode  cut-off  frequency  and:  evalua- 
tion of  f  can  be  done  by  Kuno's*  method...   The  cu.t-otff" 
frequency  is  given  by 

f  =  1/Zttt   . 

4.   Modified  Model 

The  temperature  dependence  of  diode  way  examined:  irr 
Chapter  III.   Temperature  coefficient  of  s;a.tuxa.t±an  currrejit: 
and  temperature  dependence  of  the  diode  voltage  were-  given. 
in  Eq.  3.13  and  3.14.   Equation  4.15  is  expressed  in  morre 
convenient  form  as 

CD  =  eir/m  =  (ex/mk)(I/T)   .  •  (4.16) 

For  simplification  purposes,  making  some  definitions  will  be 
useful  at  this  point  such  as  follows: 

A  -    ei/mk  4.17 

B  =  eE/mk  4..18 

C  =  Fractional  increase  of  reverse 

current  with  temperature. 

Temperature  dependent  reverse  saturation  current  increases 

exponentially  as  indicated  by  the  following  equation: 

I  T  =  I   exp[C(T  -  T  )]  4-..19 

oT    o   r  L  *•     oJ  J 


*H.  J.  Kuno,  "Analysis  and  Characterization  of  p-n  Junction 
Diode  Switching,"  IEEE  Trans,  on  Electron  Devices,  Vol..  ED- 11', 
p.  8-14,  Jan.  1964. 
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where    I   T   is   the   reverse  oirireixt  art  tre3irppena=tiirer  TT  andd  E     iss 
the   reverse   current    at   temperature  T   ..      The?  S5cr1"iit±ocnr.  ocff 
Eq.    4.19    for    C   turns    out    to   be 

C   =    ln(IoT/IQ)/CT   -    TQ)       ..  4-.2D 

For   example    1N3605   silicon  diode  foa.s^  "thj&  fixLXawijrg:  s^tes±fdra> ■ 


tions 


I      =0.05   mA   at   T     =    25°    C 
o  o 


I    T   =    50    mA   at    T   =    150°    C 


Applied  Voltage   =  —3D   volts. 


Therefore 

C   =    ln(50/0.05)/(150    -    25)    =   O..Q^/OZC     .. 

Proceeding  with  concise  form  of  Eq .  4  ..16,  the  following 
relationship  is  obtained:  * 

CD  =  A(I/T)  .  4.21 

If  we  drop  the  unity  in  the  parenthesis  in  Eq .  4.5  an  approxi 
mate  current  equation  can  be  written  as 

I  =  I  [exp(eE/mkT)]  =  I   exp(B/T)  4.22 

where  B  is  given  in  Eq.  4.18.   The  AT  change  in  temperature 
causes  the  following  modification  for  Eq.  4.22. 

IT  =  IoT  exp[B/(T  +  AT)] 

=  I   exp(C.AT)  exp[B/(T  +■  AT)] 

=  I   exp(C.AT)  exp(B/T)  exp[-E(AT)  /T(T  +-  AT)  ] 

Definition  of  another  factor  such  as 

D  =  exp(C.AT)  exp[-B(AT)/T(T  +■  AT)] 
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provides  a  very  short  expression  for  IT,  given  below: 

IT  =  D.I   .  4.23 

On  the  other  hand  AT  change  in  temperature  implies  the  fol- 
lowing modification  on  diffusion  capacitance  Cn. 

CDT  =  A.IT/(T  +  AT)  =  A.(D.I)/(T  +  AT) 

=  (A.I/T) [T.D/(T  +  AT)]=  CD  .  T.D/(T  +  AT) 

=  F.CD  4.24 

where  F  =  (T.D)/(T  +  AT).   Then  adding  and  substructing  the 
CD  to  Eq.  4.24  yields 

CDT    F,CD  +  CD  "  CD 


=  CD  +  CD  (F  -  1)  .  4.25 


In  similar  fashion 


IT  =  D.I  +  I  -  I 

=  I  +  I(D  -  1)  .  4.26 

The  temperature  treatment  of  more  sensitive  diode  parameters 
was  given  by  Eqs .  4.25  and  4.26.   The  changes  due  to  tempera- 
ture for  RR  and  R.  are  assumed  very  small  and  neglected. 
Applying  all  necessary  modifications  one  would  get  the  com- 
puter diode  model  when  temperature  vatiations  exist  as  shown 
in  Fig.  4-13. 
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'Cathode 
Figure  4-13.   Diode  Computer  Model  with  Temperature  Effects 


F.   TRANSISTOR  MODEL 


1.   General 


A  transistor  model  in  computer-aided  analysis  should 
meet  the  following  properties:   (a)  The  model  must  consist 
of  circuit  elements  that  are  recognizable  by  the  computer 
program  being  used.   (b)  The  parameters  used  in  the  model 
should  either  be  obtainable  from  published  data  or  readily 
estimated.   (c)  The  model  should  be  valid  to  high  frequencies 
as  well  as  to  low  frequencies.   (d)  The  model  should  function 
without  regard  to  the  circuit  configuration  (common  base, 
common  collector,  or  common  emitter)  in  which  the  transistor 
is  placed. 
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Note :   A  general  model  for  a>  die-trice  dices  mxt  sbDw^se  Lemhi  toe 
the  most  efficient  analysis  since  certain:  ttinjufeLesaaafire;  edieaiTejrtSi 
of  the  model  may  not  actually  be  required  by  the?  pxatrlenr  at: 
hand.   For  example,  nonlinear  or  high  frequency  elements 
may  not  affect  the  circuit  response  gx/eavtly,,  hut:  cmillcL 
easily  affect  the  solution  time  sign±ficaiT.tLy .. 

There  are  many  models  that  have  breear.  pxopxrsesdi  to: 
represent  transistors  at  low  frequencies..   Any  afT  these:  smaJ.ll-- 
signal  models  is  valid  only  when  the  transis-tor  operates? 
in  the  active  region  and  produces  negligible  ctxs tortioir 
in  the  output  signal.   Fortunately  the  autput  characteriistdcs; 
of  transistors  are  reasonably  linear  over  a  fairly  la:rge 
range,  often  allowing  output  swings  of  several  volts  with 
little  distortion.   High-frequency  models  £or  the  transistors 
are  more  complex  than  low  frequency  models,  since  additional 
reactive  elements  are  included. 

2 .   Model  without  Temperature  Effect 

The  transistor  is  a  three- terminal  device,  fabricated 
by  placing  two  p-n  junctions  in  close  proximity  to  each  other. 
The  current  through  one  forward-biased  junction  strongly 
influences  the  current  through  the  other.   The  normal  mode 
of  transistor  operation  occurs  when  the  base-emitter  junction 
is  a  forward  biased  and  base-collector  junction  is  reverse 
biased.   A  transistor  may  also  be  operated  in  the  inverse 
mode.   In  this  case,  the  collector  is  forward  biased  and  the 
emitter  is  reverse  biased.   Therefore,  transistor  can  be 
considered  as  two  junction  diodes  connected  in  series,  of. 
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which  one  is  forward  and  the  other  is  re?vre>rse?  trriasaaaii. .  TKec 
emitter  junction  and  the  collector  junction,  current,  esqnzcr- 
tions  can  be  written  as  in  single  diode  case..   By-  re^errrijrg: 
to  Section  E  and  taking  the  inverse  made  cuxreait:  gain  intn. 
account,  the  following  current  equations  may  tree  eacnrresssieadi  in. 
terms  of  diode  parameters: 

h   =  h?    +  ^VW  +  CTEfdVE/dt^  +  ^EmE/m    -  a^  V.2T 
Ic  =  Icp  ♦  (VC/RLC)  +  CTC(dVc/dt)  ♦  CDC(dVc/dt)  -  %LEF     4.28 

where  parameters  a^  and  aT  are  normal  and  inverrtexL  shtortt— 
circuit  current  gains  between  emitter  and  callercrtarr,,  capaxi:-- 
tors  CTp  and  C~p  represent  emitter-base  depletion-layer  and: 
emitter-base  diffusion  capacitances,  respectively..  WhiXe:  exp- 
and Cnr  represent  the  collector-base  depletion  and  diffusion 
capacitances,  Vp  and  V"r  are  the  emitter-base  and  collector- 
base  voltages,  respectively.   Emitter  and  collector  reverse 
saturation  currents  may  be  involved  into  the  junction  current 
equations  in  such  a  manner  given  in  the  following  equations : 

IEF  =  IES[exp(6VE/mE)  -  1]  4-.  29 

ICF  =  Ics[exp(6Vc/mc)  -  1]   .  4.30 

Ipc,  and  Irc.  are  usually  given  in  specification  sheets,  but 
the  analytical  expressions  for  Ip~  and  1^-  are  given  as 

!ES  =  WC1  "  <W  *"33: 

and 

JCS  =  W"  -    aN°l'  ^^ 
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where  Ip~  is  the  saturation  current  of  emitter  JunxrtticaT.  with 
zero  collector  current  and  I~n  is  the  s aturavtian  oianrsantt  a£ 
collector  junction  with  zero  emitter  current. 

Now  we  are  in  a  position  to  model  the  transistor 
device.   The  computer  model  uses  two  p-n  j1  unctions  mudieliexi 
like  the  diode  previously  discussed  plus  the  current  gaiir. 
factors  as  shown  in  Fig.  4-14.   The  resistances:  EL,  and  ffr 
account  for  ohmic  effects  and  can  be  approximated  as  acrnstarrrs 
in  most  instances. 

The  junction  capacitance  factors  are: 

Emitter  transition  capacitance,  CTF  =  CFnV' CvT T  ~  ^WJ  '    4^.33; 

and 

Collector  transition  capacitance,  CT~  =  C~n/(VT  -  V~)    4:„34 

where  n  =  1/2  and  n  =  1/3  for  an  abrupt  junction  and  a  linear 
graded  junction,  respectively.  Cp~  and  C^  are  determined  in 
the  same  manner  as  explained  in  previous  section  (Section  E) . 

Most  transistor  manuals  give  the  normal  and  reverse 
common  emitter  current  gains  hpp.^  and  hpp,-T>.   The  relation- 
ship between  the  four  gain  factors  is  given  by 


hFE(N)  =  V(1  "  V      aN  =  hFE(N)/(1  +  hFE(N)} 

hFE(I)  =  aI/C1  "  aI3      aI  =  hFE(I)/(1  +  hFE(I))  * 
The  diffusion  capacitances,  CDE  and  CD~  are  given  by 


4.35 


and 


CDE  =  6IEF  TE/mE  4:-36 


CDC  =  9ICF  TC/mC  4"-5T 
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where  xp  and  t„   are  emitter  and  collector  time  constants, 
respectively.   These  time  constants  can  be  evaluated  by  the 
Kuno's  Method  that  has  been  mentioned  in  Section  E. 
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Figure  4-14.   n-p-n  Transistor  Model  for  Computer  Analysis 

Most  of  the  transistor  parameters  are  supplied  by 
the  manufacturer.   However,  if  it  becomes  necessary  to  deter- 
mine the  parameters,  Refs.  18,  19,  20  and  21  are  very  useful 
for  evaluation. 

3.   Modified  Model 

The  temperature  dependence  of  transistor  parameters  was 
discussed  and  some  assumptions  were  stated  in  Chapter  III.   Addi 
tional  assumptions  are  given  below: 

(a)  aN  and  aT  are  approximately  the  same  at  temperatures 

T  and  (T  +  AT), 
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(b)  AT   <<    T, 

(c)  tf   and   x      are   the   same   at  temperr^turress  TT  ami:  (T7  +-  ATT)  ,. 

(d)  Ipp   =    Ipj,  [exp(6Vp/m.p)  ]    and  the-  same  appr odc imatiixrr 

is    for    I,,- 
CF, 

(e)  Temperature   dependence   of  EUp  ami  R>:f-  iss  ass  eaqrlaan-edi 
in   Chapter    III. 

Equation   3.17   was   written   far  cal  lector/  jumrtixnr.  ditxrfe? 
saturation   current,    and   a   similar  expression:  can-  bee  wriften 
for   emitter   jucntion   diode   saturation   current:  as 

^OT   =    ^O^1    +    AE^r^ 

=  iE0  ♦  AECArj.rEa  4i3ss 

where  I™,-  is  emitter  junction  saturation  current:  att 
fcOT  J 

temperature  (T  +  AT)  while  IFn  is  at  temperature  T..  The: 
expression  for  Ipc,  was  given  by  Eq.  4.31  for  a  specified 
temperature.   If  Ipcj  i-s  defined  as  the  value  of  Ip^  at 
temperature  (T  +  AT) ,  the  following  expression  can  be  written 

IEST  =  ^OT^1  "  aNaI^ 

=  [IEQ  +  XE(AT)  IE0]/(1  -  aNaj) 

•  hsT   =  WC1  '  Vl>  +  [AE(  T)IE0I/(1  "  aNaT) 

"  hs   +  KE  TES  4'39 

where  K~  =  A„(AT).   Equation  4.18  was  the  defining  equation 
for  the  factor  B.   A  similar  definition  can  be  used:  for 
emitter  junctions  as  follows: 

Bp  =  eVp/nipk 
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For   the   derivation   of    IEFT,   we  first  express  thee  tterrnr 
Mexp[BE/(T   +    AT)]"   in   a  more  convenient   form,  asss  fori  "Taw  ss:: 

exp[B£/(T   +    AT)]    =    exp(BE/T)    exp[-BE(AT)/r(T  -h  AX)  ] 

=    L£   cxp-CBE/T) 

where    L£   =    exp[-BE(AT)/T(T  +   AT)]. 

Thus    for    a   change   in  temperature:,.   ATI   „  the?  -FfrlTl  nwin-g 
approximate   expression   for   the    forward  emitter'  current:  is: 
obtained: 

*EFT   ^    JEST    ^*PCBE/T)ILE 

=    (IES   +    KE    1^)  [exp(BE/T)ILE 
=    IES(1    +    K£)LE   exp(BE/r) 
-    IEp    (1    +    KE)LE 

=   ME    XEF  4.40 

where   Mp   =    Lp(l    +    Kp) . 

Another  definition  given  in  Section  E  was  Eq .  4.17 
which  is  repeated  here  for  emitter  junctions  as  follows: 

A£  =  exE/MEk   . 
Then  C^p  can  be  expressed  as  CQE  =  Aplgp/T.   Therefore, 
emitter  diffusion  capacitance  at  temperature  T  +  AT  becomes 

CDET  "  AE  'EFT^1  +  AT) 

'   NE  CDE  V* 

where  N£  =  M£T/ (T  +  AT). 
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As  stated  before,  a  large  cdrangs?  irr.  hVp-,,  huweverr, ,  may* 
reflect  only  a  small  change  in  a  ami  im  nrany  ijrsftaircesf  cajr 
be  neglected  in  this  treatment.   Therefore?  the  current: 
generator  (a^I-p)  in  Fig.  4-14,  where  the  temperature  varia- 
tion is  not  considered,  can  he  written,  as?  (o^X^™^;  with. 
temperature  change  (AT) ,  so  that 

aNIEFT  =  aNME  LFiF  44.4-2- 

in  which  Eq.  4.40  has  been  applied. 

One  may  follow  the  same  derivations  given,  for  the; 
collector  junction  and  obtain  the  foLLawing;  resulltss:: 

I    =  M  I 
XCFT    C  CF 

r    =  k  c 

DCT     C   DC 

aIICFT  =  aIMC  ZCF   * 
Same  manipulation  can  be  done  for  obtaining 
temperature-caused  current  generators  and  capacitors  as 
in  the  diode  case  and  then  the  following  are  obtained: 

JEFT  =  ME  JEF  +  ZEF  "  rEF 
=  IEp  ♦  IEp  (ME  -  1) 

XCFT  =  JCF  +  XCF  ^MC  "  1-) 

r  =wr    +    r        -    r 

UDET     E   DE     DE     DE 

=  CDE  +  CDE  ^NE  "  X) 

CDCT  =  CDC  +  CDC  (NC  "  L)   *• 
Now  the  modified  computer  model  for  transistors  can  be  drawn 

as  shown  in  Fig.  4-15. 
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Figure  4-15.   n-p-n  Temperature  Effected  T 

for  Computer  Analysis 
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G.   FIELD-EFFECT  TRANSISTOR  MODEL 
1 .   General 

Although  field-effect  transistors 
used  extensively  at  present,  there  is  no  : 
describe  all  such  devices.   One  reason  is 
processes  causing  various  modes  of  operati 
another.   There  are  two  main  types,  juncti 
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1 

one 
Lated 
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gate  FET,  sometimes  called  the  "Metal-Oxide  Semiconductor" 
(MOS)  transistor  [Ref.  13].   This  would  suggest  the  use  of 
at  least  two  different  models.   A  piecewise  linear  model 
has  been  developed  by  Roberts  and  Harbourt  [Ref.  9]  that  is 
designed  for  circuit  analysis  using  the  ECAP  [Ref.  1]  pro- 
gram.  However,  piecewise  linear  models  have  some  deficiencies 
when  applied  to  general  circuit  analysis  programs  because  a 
piecewise  linear  model  can  be  made  to  approximate  the  actual 
behavior  quite  accurately  by  using  several  regions.   However, 
as  more  regions  are  added,  more  constants  must  be  provided 
and  the  awkwardness  involved  in  determining  the  correct  region 
of  operation. 

2.   Models  without  Temperature  Effects 

It  is  possible  to  make  analogy  betwpen  vacuum  tube 
(such  as  a  triode)  and  FET,  provided  that  the  cathode  is 
replaced  by  source,  the  plate  replaced  by  drain  and  grid 
replaced  by  gate.   Neglecting  the  bulk  resistances  and  taking 
the  barrier  capacitances  into  account,  the  following  equiva- 
lent circuit  can  be  obtained.   Since  the  gate  junction  is 


Gate 


GS 


♦  - 
o 


K 


*GD 


Source 
Figure  4-16. 


O 


3m  GS 


D 


Drain 


t% 


v 


DS 


-  V 
— O 


Source 
Small-Signal  JFET  Equivalent  Circuit 


53 


reverse  biased,  the  gate- source  resistance  r>^-  ami:  thee 
gate-drain  resistance  rr~  are  extremely  Large,,  and  herrcBt 
have  not  been  included  in  the  circuit  af  Fig:.  4:-Lfr..  Ileeterr-- 
mination  of  the  equivalent  circuit  parameters  will  be  dis- 
cussed later.   As  pointed  out  before  there  are?  twee  cratjegrjrie-35 
of  field-effect  transistors,  JFET  and  MDSrFET..  As?  thee  names; 
imply,  a  JFET  uses  the  characteristics  of"  a  revercssee— bdjased: 
junction  to  control  the  drain-source  current,,  while  in  the 
MOSFET  the  gate  is  a  metal  deposited  on  an  oxide  layer  and 
is  insulated  from  the  source  and  drain.   Bath  devices: 
operate  on  the  principle  of  a  "channel"'  current:  cantroTTed 
by  an  electronic  field.   The  control  mechanisms  for  the  two 
are  different,  resulting  in  considerably  different  characteris 
tics.   The  main  difference  between  the  two, is  in  the  gate 
characteristics.   The  input  of  the  JFET  behaves  like  a 
reverse  biased  diode  while  the  input  of  a  MOSFET  is  similar 
to  a  small  capacitor.   Based  on  the  above  discussion,  the 
models  shown  in  Figs.  4-17  and  4-18  are  obtained. 

In  Fig.  4-17  the  source  end  of  the  junction  is  Less 
strongly  reverse-biased  or  even  slightly  forward-biased, 
depending  on  the  instantaneous  value  of  v...   It  exhibits  the 
voltage-variable  capacitance  and  conductance  properties  of 
an  ordinary  p-n  diode.   A  lumped  model  of  this  gate-to- 
source  junction  can  be  produced  by  using  the  standard  built- 
in  model  for  a  junction  diode  to  simulate  these  effects. 
Transconductance  and  output  conductance  are  simulated  by  the 
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nonlinear  current  source,  i^c?,  which  can  be  modelled  by 
pointwise  linearization  as  described  later. 

As  stated  before,  the  gate  of  a  MOSFET  resembles 
a  parallel-plate  capacitor  with  a  thin  dielectric  layer. 
The  dielectric  results  in  a  very  low  gate  conductance,  and 
since  the  thickness  is  constant,  the  distributed  capacitance 
is  uniform  along  the  channel  and  independent  of  v^g.   This 
capacitance  can  be  split  between  the  source  and  drain  ends 
of  the  channel  and  both  lumped  capacitors  simulated  by 
standard  capacitor  model. 
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For  more  general  modelling  (including  large-signal 
behavior)  by  pointwise  linearization,  the  functional  relation- 
ship is  not  known  in  general  but  can  be  approximated  emipri- 
cally  for  a  given  device  by  measurements,  if  the  desired 
degree  of  accuracy  so  required. 


Figure  4-19.   Pinch-off 
Region  Transfer  Curve 
Showing  Pointwise 
Linearization 
Procedure 


Transconductance  g   is  the  slope  of. the  transfer 
characteristic  curve,  Fig.  4-19.   At  a  point  (V--,  ^nc)  "the 
tangent  line  is  given  by 


hs    =  Sm  (VGS  "  VGS>  +  JDS 


4.43 


where 


gm  =  (-2  InQQ/VD)[l  -  (Vrc/VD)]  . 


'm 


DSS'  P 


GS'  P- 


Inc.  in  the  pinch-off  region  is  expressed  as  in  Eq.  4.44. 


DS 


^S  =  W1  "  ^VGS/VP^ 


n 


4.44 


where  n  is  a  constant  for  a  given  unit,  whose  value  usually 
lies  between  1.7  and  2.7  (for  many  units  n  =  2  is  a  good 
assumption),  and  L^  is  the  drain  current  when  v-^  =  0. 
The  determination  of  the  necessary  model  parameters  will  be 
given  in  the  following  section. 
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3.   Determination  of  Model  Parameters 

Capacitance  determination:   For  the  model  represented 
in  Section  2,  the  common-source  short-circuit  input,  output, 
and  reverse  capacitances  are  given  in  the  following 
expressions,  respectively: 

c    =  c      +  c 
f    =  c       +  c 

OSS     DG     DS 

C     =  C 
^RSS    UDG   * 

Consequently,  the  model  capacitances  may  be  calculated  in  terms 
of  the  short-circuit  capacitances  as 

c       =  c 

UDG    ^RSS 

CGS  =  CISS  "  CRSS 

CDS  =  COSS  "  CRSS  * 
If  the  capacitances  Cy,^,  Cnc,c»  an^  CRC,o  are  no"t  supplied  by 
the  manufacturer  the  following  measurement  method  can  be  used 
to  determine  the  unknown  values.   Three  terminal  instruments 
are  the  most  useful,  since  in  addition  to  measuring  the 
capacitance  between  two  terminals,  they  are  equipped  to 
eliminate  the  effects  of  capacitance  between  these  two  termi- 
nals and  a  third  terminal.   This  is  shown  in  Fig.  4-20.   In 
Fig.  4-20  the  presence  of  C,  and  C2  does  not  affect  the 

measurement  of  C  .   In  some  instruments  the  "low"  terminal 

x 

is  grounded  and  the  terminal  connected  to  both  C,  and  C~  is 
called  the  "guard"  terminal.   Capacitance  measurement  jigs 
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are  diagramed  in  Fig.  4-21.   The  limitations  on  this 
measurement  are:   (1)  the  gate  signal  should  usually  be 
limited  to  less  than  200  mV,  (2)  frequency  should  be  kept 
between  500  KHz  and  5  MHz.   Below  500  KHz  the  chokes  must 
be  large,  and  above  5  MHz  series  resonance  and  some  other 
detrimental  effects  occur. 

Conductance  Measurement:   The  pinch-off  voltage,  V  , 
the  drain  saturation  current,  Irvcq*  the  drain  conductance, 
g~,  and  the  transconductance  g   are  four  of  the  parameters 
required  in  the  models  presented  in  this  section.   For  the 
pinch-off  region,  the  transconductance  is  given  by  Eq .  4.45. 


gm  ■  c^WVr1  -   (VGS/Vp)] 


4.45 


if  V^r,  =  0  then  g   is  maximum, 
GS  &m  ' 


g 


max 


or 


2  ^SS^P 


:/g, 


VP    ~2  XDS 


The  experimental  circuit  for  determination  of  g     is  given 

r  "max    ° 

in  Fig.  4-22. 

The  determination  of  gn  can  be  done  by  a  circuit 

as  shown  in  Fig.  4-23,  in  which  g„ (micromho)  =  100  V-CmV) 

when  V,  =  1.0  volt. 
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Figure  4-20.   Three  Terminal  Measurement 
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Figure  4-21.   (a)  CT~S  Measurement,  (b)  Coss  Measurement, 

(c)  Cpoq  Measurement  Jigs 
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Current  Measurement:  * 

Inc;c,  is  the  easiest  determined  parameter  as  far  as  the 
measuring  circuit  is  concerned.   The  following  configuration 
can  be  set  for  determination  of  I^cc 
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Figure  4-24.   Circuit  for  Measurement  of  I^cc 
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4.   Modified  Model 

As   pointed   out    in   Chapter   III,,   tire?  nrasjcrrr  efEfecrss  off 
temperature    on   a   FET   are    to   alter  tlie   average:  cframrerll  crjrrduc:-- 
tivity,    a,    and    the    gate- channel    junction  barrier  pcrterrtiaL,    \\j^. 
Since   a   decreases   with    increasing   tempera ture,.   the?  effect:  of; 
this    to   decrease    the    drain   current   as   "the  terarfenratiire:  risjes;. . 
On   the    other   hand,    [  t[>  [    decreases  with,   imcreasing:  temperatxire^ . 
causing    the   magnitude   of   the   total   gate-source  pertenTdLall  tro:  bee 
reduced,    hence    increasing    the   drain   current.. 

Next,    consider   the    three   terminal  FET  tn  he  operated-  in- 
the   saturation   region  with    a   constant   gate-s^crurcs  emf:..     There?-- 
fore,    Inc.    =    1^  (ij/,a)  .      The   differential   equation:  fxrr  the   drain', 
current    is    given    [Ref.    8]    by   the    following:  relationship::: 

dIDS/dT    =    (3IDS/3ifO(dU;/dT)    +    (3IDS/3a)(dcx/dT) 

=  gm  (dip/dT)  +  (IDS/a)(da/dT)  4.46 

since  31^-/3^  =  g   and  In~  is  proportional  to  a  at  a  given 
temperature  T.   A  reasonable  approximation  to  the  mobility- 
temperature  function  [Ref.  2]  is 

M  =  P0(T/To)-n 
where  n  depends  on  the  impurity  concentration.   Therefore,  after 
necessary  manipulation  the  following  equation  is  obtained: 

(1/cr)  (da/dT)  =  -n/T  4.47 

If  AT  <<  T  then  d^/dT  can  be  approximated  by  ±2.2  mV/°C 
(depending  upon  whether  the  transistor  is  p-channel  or 
n-channel,  respectively)  at  low  source  current  level!..  For:  many 
purposes  n  may  be  approximated  to  2  (actually  it  ranges  from  1.7  to 
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2.7).      Taking    the   above   app tow.t matx cure  rmtDx  adsjzmurtt  ttiee 
Eq.    4.48    can   be   obtained  far  drain,  Gmrren.tL  charrgjee  w.itir. 
respect    to    temperature. 

dIDS/dT   =    2.2(mV/°C)    gm    -  n    L^/T  4-. 4-8 

It    is    easily   seen   from  Eq.    4.45  tlxa7t  the;  czmdditiiorr.  for:  zetd: 
temperature    coefficient    is    ap pr ox±mavte:Ly •  giverr  by 

IDS/gm   =    Z.Z(T/n)    mV    ..  4445? 

When    the    change    in   temperature    is  very   smaXX  compared:  to: 

the    given    temperature,    T    (AT   <<    T)  „   Eq..   4^..48E  can.  bee  expTessed 

in   the    following    form: 

AIDS/AT    S    2.2(mV/°C)    gm    -   n    L^/T 
then 

AIDS    a    2.2    gm    (AT)    -    (n/T)    IDS(AT)     .  4,5D 

4 

Therefore,  one  may  represent  Eq.  4.50  using  two  current 
generators  with  current  directions  opposing  to  each  other. 
Applying  the  Eq.  4.50  to  Fig.  4-17  the  FET  computer  model 
is  obtained  for  a  given  temperature  T  and  the  temperature 
change  AT  <<  T,  as  shown  in  Fig.  4-25.   Tf  AT  is  large,  one 
has  to  make  iterative  computations  by  taking  proper  AT 
increments . 

The  necessary  model  parameters  may  be  determined  as 
explained  in  the  previous  section.   The  currents  I,  and  I~ 
caused  by  the  temperature  can  be  expressed  in  terms  of  known 
parameters  as 

Ix  =  2.2  gm  (AT)  4^51". 
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2  =  (n/T)  IDS  (AT) 


4.52 


To  be  dimensionally  correct  g   should  be  in  millimhos, 

°m  ' 

temperature  T  and  (AT)  in  degrees  Kelvin. 
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Figure  4-25.   Lumped  JFET  Comrmter  Model  with  Temperature 

Effect. 

The  enhancement- type  MOSFET  has  a  characteristic 
threshold  voltage,  V"T,  instead  of  the  pinch-off  voltage,  Vp, 
but  the  equations  are  of  the  same  general  form:   for  0  <  v 


-  VT  <  Vnq  (pinch-off  region)  inQ  =  K(v 


DS 


DS 


GS 


GS 
2 
VT)  .   Here  K  is 


a  constant  and  is  a  function  of  channel  geometry,  doping,  and 

effective  carrier  mobility.   In  this  section  an  analysis 

is  given  for  the  temperature  dependence  of  vrc.  when  the 

current 'is  maintained  constant  and  saturation  conditions 

are  present.   It  is  assumed  that  interface  charge  Q~~  is 

independent  of  temperature  and  that  the  effects  of  a  finite 

drain  resistance  can  be  neglected.   Further,  it  is  assumed 

that  the  only  temperature  dependent  quantities  are  v~<,  (DC 

gate-source  voltage)  ,  <£„  (Fermi  potential)  and  y   (effective 

r  n 
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surface  mobility) .   The  temperature  coefficient  as  a 
function  of  v^  is  given  [Ref.  8]  by 

dvGS/dt  =  (iD/gm)(-l/yn)(dyn/dT) 

4.53 
-  (d<f>F/dT){[-2(vGS  -  VT)/vDS]  +  1} 

Since  the  first  and  the  second  terms  are  both  positive, 
it  is  evident  that  there  may  exist  a  realizable  zero  tempera- 
ture coefficient  point.   Supposing  that  the  mobility  is  pro- 
portional to  T   *   and  letting  n  =  2,  one  can  obtain  (by 
taking  certain  known  factors  and  coefficients  into  account) 
the  following  expression  for  temperature  coefficient  [Ref.  8] . 

dvGS/dt  =  (3.0/T)(vGS  -  VT) 
-  [(0.605  -  <f>F)/T]{[2(vGS  -  VT)/vDS]  -1}  . 

As  pointed  out  by  Giralt  et  al.  [Ref.  23] , 'a  good  approxima- 
tion to  the  behavior  is  given  by 

d|vGS|/dT  =  h  -  2.5|vGS  -  VT|mV/°C   ,         4.55 

where  VT  is  the  measured  threshold  voltage  obtained  from  the 
straight  line  extrapolation  of  an  I^/g   plot  on  the  v^-  axis, 
and  h  is  a  constant  whose  value  is  somewhat  dependent  on  the 
subsrate  doping  but  normally  lies  between  2  and  6. 

If  AT  <<  T  then  the  Eq.  4.55  can  be  written  as  in 
Eq.  4.56. 

A | vGS |  =  h(AT)  -  2.5|vGS  -  VT|(AT)mV  .       4.56 

On  the  other  hand  when  temperature  changes  as  much  as  AT 
Eq.  4.43  becomes 
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=     gm[(V,Q     +      |Avpe|)      -     VfiR]      +     I 


DST        6mLV"GS 
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DS 


=    2      v^c    +    g      Av~c 
6m      GS        6m'       GS 
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6m      GS 


DS 


■    inc    +    gjAv 


DS 


m 


GS 


4.57 


This  change  Avrc,  can  be  obtained  either  from  Eq.  4.54  or 
Eq.  4.56  depending  upon  the  given  (or  known)  information. 

The  term  (g  |Av~0|)  can  be  defined  as  temperature  resultant 

V6m '   GS ' J  r 

current  generator  than  temperature  dependent  computer  model 
for  MOSFET  obtained  as  shown  in  Fig.  4-26. 
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Figure  4-26.  Lumped  MOSFET  Computer  Model  with  Temperature  Effect 
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V.   CONCLUSION'S 

A.   GENERAL 

It  has  been  shown  that  temperature  variations;  have 
various  effects  on  different  electronic  components..  Resis- 
tance of  a  resistor  increases  or  decreases  with,  temper.atiire. 
depending  on  the  material  used.   Inductance  off  air.  inane-Tor: 
may  or  may  not  change  with  temperature  depending:  on  its;  physi;-- 
cal  dimensions  and  the  material  used-   Capacitance  of:  a: 
capacitor  changes  slightly  with  temperature..   However,.  this; 
small  change  can  be  extremely  important  in.  same  situations; 
(e.g.,  as  in  resonance  circuits).   Forward  diode  current's-  and; 
diffusion  capacitance  increase  with  increasing  temperature. 

4 

The  effect  on  a  transistor  is  quite  complex,  and  some  simpli- 
fying assumptions  are  listed  in  Chapter  IV.   The  junction 
forward  current  and  diffusion  capacitance  increased  with 
increasing  temperature.   Reverse  saturation  current  also 
increased.   Two  types  of  field-effect  transistors  were 
examined  separately.   It  was  found  that  with  the  proper 
operating  point,  the  temperature  effect  could  be  eliminated. 
If  this  operating  point  cannot  be  attained,  the  drain  current 
changes.   Therefore,  the  effects  discussed  above  should  be 
applied  to  the  computer  models  when  one  attempts  to  analyze 
a  circuit  which  is  composed  of  electronic  components.. 

The  models  proposed  here  allow  simulation  olr  the  behavior: 
of  several  circuit  components.   For  a  reasonably  accurate 
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simulation,  the  parameters  mssqiiitiired:  finr  titer  moiieds?  wearee  few 
easily  determined.   The  actual:  ciircaiit:  irr.  which:  tire:  component: 
is  to  be  used  is  a  major  factor  in.  thee  cirorice:  off  as  model  I  and: 
the  use  of  these  models  must  be  done  with  an  understanding  of 
the  assumed  operating  regions  ami  cdiaractjea^isitirsf  of:  thet 
components  themselves.   Thie?  panamete3!S5  of"  titec  modeil  must"  also: 
be  chosen  to  suit  the  particular  appliication , .  because:  evenr. 
the  best  models  can  only  approximate  the  behavior:  of:  th.e: 
actual  devices. 

Emphasis  was  placed  on:  empirical  deie-rminatii)n  off  parameters ; 
however,  these  parameters  may  also:  he  datea^nrined:  fiomr  thev 
material  characteristics  and  structural  geometry  of:"  the.  devices 
being  modelled. 

These  models  are  useful  in  analysis  and  design  of  elec- 
tronic circuits.   Also,  since  the  model  parameters  may  be 
calculated  from  information  on  the  construction  of  the  device, 
it  should  be  possible  to  predict  the  behavior  of  a  type  of 
component  before  it  is  fabricated. 

The  models  are  general  so  that  they  are  compatible  with 
most  of  the  electronic  circuit  analysis  programs  in  use  today. 

B.   SUGGESTIONS  FOR  FURTHER  STUDY 

1 .   Modification  of  Models  with  High  Frequency  Effects 
High-frequency  models  for  the  components  are  more 
complex  than  low  frequency  models,,  since  additional  reactive 
elements  must  be  included..   As  a  suggestion:  for:  further  study ,  . 
one  might  consider  the  effects  of  high  frequency  effects 
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upon   these   components   under  normal!  ajBea^tdzrgr  temperatures 5 
and   temperature   as   a  variafrLe.. 

2.  FET   in   Triode    Reg-ion 

As    far   as    the    FET   is   concerned,,   the   developments 
were   done    for   the  saturation,  region,  oqierration . .     Th.ee  other: 
very   important   application  is   the  trriaxfe-re^irmr.  operation:. . 
In   triode    region,    {G    <  ^rjcr  <:  v'c<z  ~  ^pr)    tfree  g™--  is?  expressed 
as 

sDS  =  eiDSS/vj)[z(vGS  -  vF)  -  zv^  t  531: 

At  an  arbitrary  paint  (v~nc-„  O  arr.  the:  output:  characteristic: 
curves,  the  pointwise  linear  approximation  to:  the:  output: curve 
is  the  straight  line  i^   =  gDS  ^ns  -  VQS)  <-  1^  where,  gDS; 

is  the  slope  at  CVnc»  *n<^  g"iven  with  Eq.  55.JL..  Then  it" 

« 
should  be  simple  to  develop  models  for  triode  region  for 

various  FET's  by  taking  the  temperature  effects  into  account. 

3 .  Some  Other  Device  Models 

One  might  also  develop  computer  models  for  additional 
devices  depending  upon  the  methods  stated  in  this  paper  with 
temperature  or  high  frequency  effects  in  mind. 

(a)  Zener/avalanche  diode  model:  Since  the  avalanche 
diode  is  nearly  identical  to  an  ordinary  diode  for 
operating  voltages  greater  than  VR  (breakdown  voltage) , 
it  should  be  possible  to  modify  the  Fig.  4-9,  diode 
computer  model  to  include  the  avalanche  phenomenon. 

(b)  Tunnel  diode  model::   The  unique  properties:?  of:  tunnel  1 
diodes  which  are  most  often  utilized  are  its  negative 
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resistance  and  its  very  rapid  switching  from  Vp 

(peak  voltage)  to  a  point  in  the  forward  conduction 

region.   Most  tunnel  diode  circuits  are  designed  around  the 

peak  point  and  valley  point.   Any  tunnel  diode  model  should 

accurately  approximate  these  points  as  well  as  modelling 

other  important  properties.   It  is  relatively  easy  to 

measure  the  peak  point,  valley  point  and  the  maximum 

slope  in  the  negative  conductance  region.   There  is 

very  little  resistance  to  current  flow  when  reverse  bias 

is  applied  to  the  junctions;  i.e.,  with  reverse  bias, 

the  tunnel  diode  is  nearly  a  short  circuit. 

(c)  Unijunction  transistor:   On  the  emitter  characteristics 

the  peak  point  is  defined 

by  the  expression 

V      =   v   +    aVnil   where   v 
p  BB 


is  defined  by  the  IR?  =  0 
curve.   The  quantity  a, 
called  the  intrinsic 
standoff  ratio.   When  Ip  becomes  greater  than  I   (peak  point 
current) ,  the  process  of  conductivity  modulation  occurs  in 
the  semiconductor  material  between  emitter  and  base-1, 
causing  the  negative  resistance  characteristic.   As  I£  is 
increased  well  beyond  I   there  is  also  a  slight  conductivity 
modulation  of  the  B,  -  B~  region  resulting  in  an  increased 
IR2.   A  model  which  will  adequately  describe  the  static 
behavior  of  a  unijunction  transistor  can  be  formulated. 
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(3)    Silicon   controlled   rectifiers::      Since:  SFrTiiionir,  CTcntrrxL 
rectifier    (SCR)    is    a   four-layer T   tn:ree-jumrt±ajr.  dev.ixec  w.ithr 
only   three    leads,    it    is    difficult    to   describe  the  internal! 
behavior.      The    two- trans  is  tor   analogy   praviides-  3l  g.aoxL 
intuitive    feel    for   the   operation,,  hut:   it  iss  dazEExciilJt:  ti>: 
obtain   the   transistor  parameters   since  frase:— II  ami  axLLe  etc  r?  2! 
are   not    accessible.      It   should  be  possitrle:  to:  represent: 
an   SCR  by   empirically   relating   the    internal  behavior  to: 
terminal    characteristics    even   the   physical!  processes-  are 
somewhat   more    complex   than   for  same   of  the  devices:  dis.-- 
cussed   previously. 
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